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Abstract 

The bottom quark IS' mass, M^^ , is determined using sum rules which relate the masses and 
the electronic decay widths of the T mesons to moments of the vacuum polarization function. 
The IS* mass is defined as half the perturbative mass of a fictitious ^Si bottom-antibottom quark 
bound state, and is free of the ambiguity of order Aqcd which plagues the pole mass definition. 
Compared to an earlier analysis by the same author, which had been carried out in the pole mass 
scheme, the IS* mass scheme leads to a much better behaved perturbative series of the moments, 
smaller uncertainties in the mass extraction and to a reduced correlation of the mass and the strong 
couphng. We arrive at M^^ = 4.71 ± 0.03 GeV taking as(Mz) = 0.118 ± 0.004 as an input. From 
that we determine the MS mass as rfibijnb) = 4.20±0.06 GeV. The error in mb(jnb) can be reduced 
if the three-loop corrections to the relation of pole and MS mass are known and if the error in the 
strong coupling is decreased. 



PACS numbers: 14.65.Fy, 13.20.Gd, 13.20.Gv. 



1 Introduction 



The determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements is one of the main 
goals of future B physics experiments. A sufficiently accurate determination of the size of the CKM 
matrix elements and their relative phases will lead to a better understanding of the origin of CP 
violation, the structure of the weak interaction, and, possibly, to the establishment of physics beyond 
the Standard Model. For the extraction of the CKM matrix elements from inclusive B decay rates, 
particularly of Vcb from the semileptonic decays into D mesons, an accurate and precise knowledge 
of the bottom quark mass is desirable due to the strong dependence of the total decay rate on the 
bottom quark mass parameter. 

Non-relativistic sum rules for the masses and electronic decay widths of T mesons, bottom- 
antibottom quark bound states having photonic quantum numbers, are an ideal tool to determine 
the bottom quark mass: using causality and global duality arguments one can relate integrals over 
the total cross section for the production of hadrons containing a bottom-antibottom quark pair in 
e~^e~ collisions to derivatives of the vacuum polarization function of bottom quark currents at zero- 
momentum transfer. For a particular range of numbers of derivatives the moments are saturated by the 
experimental data on the T mesons and, at the same time, can be calculated reliably using perturbative 
QCD in the non-relativistic expansion. Because the moments have, for dimensional reasons, a strong 
dependence on the bottom quark mass, these sum rules can be used to determine the bottom quark 
mass to high precision. Within the last few years there have been several analyses at NLO and NNLO 
in the non-relativistic expansion, where it has been attempted to extract the bottom quark pole mass 
from these sum rules, fl], I, I, I, IQ The results for the pole mass given in these analyses vary rather 
strongly in their central values and their precision. This reflects the fact that the concept of the pole 
mass is ambiguous to an amount of order the typical hadronization scale Aqcd 0) §]• From a technical 
point of view the large variations in the extracted pole mass values arise from large correlations to the 
value of the strong coupling a^, the theoretical uncertainties coming from large NNLO corrections to 
the moments, a strong residual dependence on the renormalization scale, and a different interpretation 
and treatment of these sources of uncertainties in the various analyses. From a conceptual point of 
view at least part of the former issues can be traced back to the fact that in the pole mass scheme the 
theoretical expressions for the moments are sensitive to scales which are below the scales characteristic 
to the non-relativistic bottom-antibottom quark dynamics encoded in the moments. This sensitivity 
to low scales, which becomes stronger at higher orders of perturbation theory and which is called the 
"renormalon" problem of the pole mass definition, is an artifact of the pole mass scheme, and does 
not exist if a mass definition is employed which is not ambiguous to an amount of order Aqcd |9|, P^ . 
Such mass definitions are referred to as "short-distance" masses. (Short-distance masses can still have 
ambiguities of order Aq^^j-, divided by the heavy quark mass.) In general, these mass definitions have 
a nicely behaved perturbative relation to the MS mass. This means that the perturbative relation 
between the short-distance masses and the MS mass can be expected to be convergent for all practical 
purposes (i.e. as long as we only deal with a few low orders of perturbation theory).]^ However, 

^ In Ref. a sum rule analysis has been performed which is not compatible with the non-relativistic velocity counting 
rules. 

^ The words like "convergent" and "nicely behaved" , which are used in this work to describe the perturbative relations 
between short-distance masses, do not have a strict mathematical meaning. They are merely chosen to distinguish from 
the situation one finds if the pole mass is related to a short-distance mass. We would describe the latter situation with 



2 



not any short-distance mass is equally well adapted to improve the situation for the T sum rules 
because the correlations to the strong coupling and the strong dependence on the renormalization 
scale observed in the pole mass scheme already exist in the leading order theoretical expressions for 
the moments and are therefore not related to the problem of an increased infrared sensitivity. Thus, it 
is advantageous to use a specialised short-distance mass for the mass extraction from the T sum rules, 
which eliminates, as much as possible, correlations and dependences to other parameters and which 
stabilises the perturbative expansion of the moments. Such a specialised short-distance mass can be 
either used as a mass definition in its own right, or, in a second step, related to other short-distance 
masses like the MS mass, which can be regarded as a specialised short-distance mass designed for high 
energy processes. In this step a part of the correlations to the strong coupling are expected to come 
back. However, correlations to other parameters like renormalization scales are eliminated, which can 
lead to a reduction of uncertainties. In addition, the perturbative relations between short-distance 
masses is well convergent and possible sources of correlations which arise in these relations can be 
easily identified. 

In this paper we extract the 15 bottom quark mass from the T sum rules based on the 
theoretical expressions of the moments which we have calculated analytically at next-to-next-to-leading 
order (NNLO) in the pole mass scheme in a previous publication Q. In the IS" scheme the bottom 
quark mass is defined as half the perturbative mass of a colour singlet bottom-antibottom quark 
J^^ = 1~, ^Si bound state. The 15 mass is a short-distance mass and has been used previously 
to parameterise inclusive B mesons decays leading to a considerable reduction of the perturbative 



corrections. 12 1 In the 15 mass scheme we find a reduction in the size of the large NNLO corrections 



to the theoretical moments observed in the pole mass scheme, and a much weaker dependence of the 
moments on ag and the renormalization scale governing the non-relativistic bottom-antibottom quark 
dynamics. This results in a weak correlation of the 15 mass to the value of the strong coupling 
and to much smaller uncertainties in the mass determination compared to our earlier analysis in 
the pole scheme {M^^ = 4.71 ± 0.03 GeV in this work versus Mj^"^'' = 4.9 ± 0.1 GeV in Ref. |] if 
as{Mz) = 0.118 lb 0.004 is taken as input) using exactly the same statistical analysis and the same 
way to treat theoretical uncertainties. 

By construction, the 15 mass is half the perturbative contribution of the T(15) meson mass 
{Mr (IS) = 9460.37 ± 0.21 MeV H). In Refs. 0, this relation has been used explicitly, and an 
uncertainty of ±50 MeV has been assigned to the value of M^^ based on conservative general arguments 
to estimate the size of non-perturbative effects in the T(15) mass. In this work this relation is not 
used and M^^ is treated as a fictitious mass parameter which is determined solely from the sum rule 
analysis. The close proximity of M^^ determined in this work and My(i5)/2 = 4.730 GeV is non- 
trivial because higher radial T meson excitations and the bb continuum have a significant contribution 
to the sum rule. Thus the analysis in this work provides an independent quantitative cross check of 
the arguments about the size of non-perturbative effects in the T(15) mass made in Refs. |11, 12|. 

The program of this paper is as follows: in Sec. |2| we introduce the 15 mass. In Sec. ^ we 
briefly review the T sum rules and discuss how the analytic results obtained in Ref. |Q] in the pole 
scheme are modified in the 15 mass scheme. In Sec. ^ the statistical analysis is explained and the 
results of the determination of the 15 bottom quark mass are shown. In Sec. ^ we discuss the upsilon 

the words "not convergent" and "badly behaved". From the mathematical point of view, of course, all series shown in 
this work are asymptotic. 
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expansion [11, 12| which is needed to relate the IS* mass to other short-distance mass definitions, and 
we determine the MS mass. In Sec. ^ we determine the value of other recently proposed specialised 
short-distance masses which can be used in the T sum rules, and Sec. |^ contains the conclusions. 



2 The IS Bottom Quark Mass 

The IS" bottom quark mass is defined as half the perturbative mass of a J^*^ = 1 , ^Si bottomonium 
ground state. Expressed in terms of the pole mass the IS" mass at NNLO in the non-relativistic 
expansion reads (a^ = as{fi)), |14, p[ 
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The constants /3o and /3i are the one- and two-loop coefficients of the QCD beta function and the 
constants ai [|l^, |l6| and 02 [|T7| , [T8|| the non-logarithmic one- and two-loop corrections to the static 
colour-singlet heavy quark potential in the pole mass scheme, yC"""'.^ The charm quark is treated as 

^ The constant 02 was first calculated in Ref. Recently an error in the coefficient of the term oc ii'^C\ was 

corrected in Ref. [ |l7[ |. Although this leads to a reduction of the value of a2 by a factor of two, the change turns out to be 
irrelevant for the mass determination in this work, because the influence of 02 is, by construction, strongly suppressed 
in the IS scheme. 
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massless. In Eq. (|l|) we have labelled the contributions at LO, NLO and NNLO in the non-relativistic 
expansion by powers e, and e'^, respectively, of the auxiliary parameter e = 1. The expansion in 



terms of the parameter e is called the upsilon expansion |11, |l2| and will be relevant when the IS is 
related to the MS mass. We will come back to this issue in Sec. ^ 

In the framework of the non-relativistic power counting scheme A^O) ^nlo and An^iq are 
of order v"^, and respectively. In order to implement the IS" mass into the analytic expressions for 
the moments in the pole mass scheme we have to invert relation in the non-relativistic framework, 
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I 1 + A^°{as) + A 
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A™0(M,i^a„/.) + (A^«(a, 



(7) 



The terms in the brackets on the RHS of Eq. (^) represent the NNLO contributions. We would like to 
emphasise again that M^^ is a short-distance mass because it does not suffer from the ambiguity of 
order Aqcd like the pole mass. This is because the IS mass contains, by construction, the half of the 
total static energy can be proven to be ambiguity-free at order Aqcd- ||9|, |1^ 



3 The T Sum Rules 

In this section we briefly review the basic concepts involved in the T meson sum rules. We outline 
the calculation of the moments in the pole mass scheme which has been carried out in Ref. |^, and 
we describe how the moments have to be modified in the IS scheme. We will not present any details 
on the original computations carried out in the pole mass scheme and refer the interested reader to 
Ref. i. 

The sum rules for the T mesons start from the correlator of two electromagnetic currents of 
bottom quarks at momentum transfer q, 

n^.(9) = -i y'dxe*''-(0|rj^(x)j^(0)|0), (8) 

where 

=6(x)7^6(x), (9) 

and the symbol b denotes the bottom quark Dirac field. The n-th moment P„ of the vacuum polar- 
ization function is defined as 



(10) 

g2=0 



where Qh = —1/3 is the electric charge of the bottom quark. Due to causality the n-th moment P„ 
can also be written as a dispersion integral 

no 

ds „, . 
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where 

R{s) = ^ (12) 

is the total photon mediated cross section of bottom quark-antiquark production in e~^e~ annihilation 
normalised to the point cross section apt = 47rQ;^/3s, and s the square of the centre-of-mass energy. 
The lower limit of the integration in Eq. (lTl| ) is set by the mass of the lowest lying resonance. Assuming 
global duality the moments Pn can be either calculated from experimental data on R or theoretically 
using perturbative QCD. 

The experimental moments R^^ are determined using latest data on the T meson masses, 
^T{n5)) ^-nd electronic decay widths, T-f^^nS)^ for n = 1, . . . , 6. The formula for the experimental 
moments used in this work reads 

g oo 

^ /^ww. (13) 



-em j^^i ^''-'T(fc) 



BB 



and is based on the narrow width approximation for the known T resonances. Oem is the electro- 
magnetic coupling at the scale 10 GeV. Because the difference in the electromagnetic coupling for the 
different T masses is negligible we chose 10 GeV as the scale of the electromagnetic coupling for all 
resonances. The continuum cross section above the BB threshold is approximated by the constant 
Vc = 1/3, which is equal to the born cross section for s oo, assuming a 50% uncertainty^ 

W(s) = r, (1 ± 0.5) . (14) 

For n > 4 the continuum contribution is already sufficiently suppressed that a more detailed description 
is not needed. For a compilation of all experimental numbers used in this work see Tab. 

A reliable computation of the theoretically moments based on perturbative QCD is only pos- 
sible if the effective energy range contributing to the integration in Eq. (0) is sufficiently larger than 
^QCD ~ 0(200 — 300 MeV) (2^]. For large values of n one can show that the size of this energy range 
is of order Mfj/n. This means that n should be chosen sufficiently smaller than 15 — 20. To suppress 
systematic theoretical uncertainties as much as possible we take Umax = 10 as the maximal allowed 
value for n. However, it is also desirable to choose n as large as possible in order to suppress the 
contribution from the 66 continuum to R{s) above the BB threshold, which is rather poorly known 
experimentally. In other words, one has to choose n large enough that the bottom-antibottom quark 
dynamics encoded in the moments P„ is non-relativistic. Because the effective size of the energy range 
contributing to the n-th moment is of order Mfj/n, the mean centre-of-mass velocity of the bottom 
quarks in the n-th moment is of order 

w ~ tts ~ , (15) 



where f ~ is characteristic for perturbative non-relativistic quark-antiquark systems. In our anal- 
ysis we choose nmm = 4 as the minimal value of n to ensure the dominance of the non-relativistic 
dynamics in the moments. For the values of n employed in our analysis the non-perturbative contri- 
butions coming from the gluon condensate are at the per-mill level and negligible. 

* We take the opportunity to point out a typo in Eq. (77) of Ref. where a factor 1/3 is missing. This typo only 
exists in the text of Ref. ^ and is not contained in the numerical codes. 
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nS 


IS 


25 


35 


Mns/[GeY] 


9.460 


10.023 


10.355 


rn5/[keV] 


1.32 ±0.04 ±0.03 


0.52 ±0.03 ±0.01 


0.48 ±0.03 ±0.03 


nS 


4S 


55 


65 


M„5/[GeV] 


10.58 


10.87 


11.02 


r„5/[keV] 


0.25 ±0.03 ±0.01 


0.31 ±0.05 ±0.07 


0.13±0.03± 0.03 


d-^ = Q-^(10 GeV) = 131.8(1 ± 0.005) , (^)bb = 2 x 5.279 GeV 



Table 1: The experimental numbers for the T masses and electronic decay widths used for the 
calculation of the experimental moments P^^- For the widths the first error is statistical and the 



second systematical. The errors for the partial widths of T15 and T25' are taken from Ref. |19]. All 



the other errors are estimated from the numbers presented in Ref. [13|. The small errors in the T 
masses and the BB threshold {y/s)BB ^^^^ neglected. 



In Ref. we have used non-relativistic QCD (NRQCD) as formulated in Ref. [21| to determine 
the theoretical moments in the pole mass scheme at NNLO in the non-relativistic expansion, 
which included all corrections up to order 1/n, asj^fn and o?^ with respect to the expressions which 
are leading order in the non-relativistic expansion. Taking into account the power counting shown in 
Eq. ([1^) the dispersion integration for the theoretical moments Pl^ at NNLO takes the form 

pth ^ 1 /■ ^-^ [ ^ n\{\ ^ I ^ R'''' (E) (16) 



bind 



where E = y/s — 2M^°^'^ and -Ebind is the binding energy of the lowest lying resonance, i.e. -Ebind = 
2{Ml^-Ml°^''). The exponential form of the LO non-relativistic contribution to the energy integration 
has to be chosen because E scales like v"^ ~ 1/n. The result for the theoretical moments at NNLO 
can be cast into the form 

pth _ ^NcQl a/tF J / ^hard /^fac / n\ / /^soft /^fac / n\ 

+ 0n,2{as{fisoit)) I , (17) 

where Qn,i describes the contribution to the moments coming from the dominant non-relativistic 
current correlator involving two dimension three ^5i NRQCD currents, and Qn,2 the contribution 
coming from the NNLO current correlator involving one dimension three and one dimension five 
^5i NRQCD current. Ci contains the short-distance corrections to the dimension three currents 
up to order O(a^). The corresponding short-distance correction to Qn,2 is not needed because the 
latter is already of NNLO. In Eq. ( |T7| ) we have indicated the dependence of the moments on the 
various renormalization scales used in Ref. . /igoft is the renormalization scale of the strong coupling 
governing the non-relativistic dynamics and //hard the scale of the strong coupling in the short-distance 
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Moment 


Mi^(Mf' ")/[Gey] 


as{Mz) 




4.6 


4.8 


5.0 


0.11 


0.12 


0.13 


Pf /[10-8GeV-8] 


0.27(0.39) 


0.19(0.28) 


0.14(0.20) 


0.17(0.22) 


0.20(0.30) 


0.26(0.46) 


Pf /[10-12 


0.28(0.47) 


0.17(0.29) 


0.11(0.18) 


0.14(0.21) 


0.18(0.31) 


0.24(0.54) 


Pg / [iU GeV J 


0.33(0.64) 


0.17(0.33) 


0.09(0.17) 


0.14(0.22) 


0.18(0.36) 


0.25(0.70) 


oth /r-i n-20 noAf-20i 
Plo/[iU GeV J 


0.41(0.90) 


0.18(0.39) 


0.08(0.17) 


0.14(0.25) 


0.19(0.44) 


0.27(0.94) 




1.64(6.85) 


0.30(1.26) 


0.06(0.25) 


0.22(0.63) 


0.33(1.52) 


0.53(4.72) 




as{Mz) = 0.118 


Mfti^ = 4.8 GeV 
(M^'"'" = 4.8 GeV) 




;Usoft = 2.5 GeV , /ihard = Aifac = 5 GeV 



Table 2: The theoretical moments F^^ at NNLO in the IS (pole) mass scheme for n = 4, 6, 8, 10, 20 
and fixed figoit = 2.5 GeV and ^hard = Wac = 5 GeV for various values of M^^ (M^°'^) and as(M^). 
The two-loop running has been employed for the strong coupling. The values for P^^ in the pole mass 
scheme are slightly different from the numbers shown in Ref. [^ because of the correction in a2- 



coefficient Ci. The scale //fac is the factorisation scale which separates non-relativistic and short- 
distance momenta. The three scales arise in the calculation of Ref. [Q| as a consequence of the use of 
a cutoff-like regularization for the UV divergences which arise in NRQCD Feynman diagrams, while 
the strong coupling is renormalised in the MS scheme. Formally the moments are invariant under 
changes of these scales at NNLO. Prom the numerical point of view, however, they are not, because 
the dependence on the scales, in particular the soft scale, only cancels partially at finite order of 
perturbation theory. In the statistical analysis for the bottom mass determination all three scales 
are varied independently in order to estimate theoretical uncertainties. Eq. (|T7| ) also displays the 
dependence on the pole mass. The short-distance factor Ci and Qn,i only depend on the pole mass 
through the logarithm of the ratios of the renormalization scales and the pole mass, which originate 
from the running of the strong coupling and the NRQCD UV divergences. The latter dependences 
only arise at NLO and NNLO and do not lead to any modifications in the 15" scheme because the 
difference between 15" and pole mass is of order ~ ~ 1/n in the non-relativistic expansion. The 
most important pole mass dependence is the overall factor (M^"''^)-^". In the IS* scheme this factor 
reads 

1 1 



+ n 



exp ( - 2n A^°(a,(^,oft))) 1 1 - 2 n A^^°(Mbi^, a,(^soft), /^soft) 

2 „ 1 o / „, „ 1 a N 2 



(A^°(a,(/Xsoft))) - 2A^^^°(M,i^, a.(;Usoft), //soft) + 2n(A^^°(Mi^, a.(^soft), /isoft) 

(18) 

at NNLO in the non-relativistic expansion using the relation given in Eq. (^). Because the IS mass 



is defined purely from the non-relativistic dynamics, the soft scale has to be employed in Eq. (18). 



Thus, the transition to the IS" mass scheme leads to a simple rescaling of the theoretical moments. 
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Moment 


/igoft / 


GeVl 


/"hard /[GeV] 


/^fac/ 


GeVl 

J 




1.5 


3.5 


2.5 


10.0 


2.5 


10.0 


Pf /[lO-^GeV-S] 


0.33(0.51) 


0.17(0.23) 


0.17(0.25) 


0.21(0.30) 


0.21(0.31) 


0.16(0.22) 


Pi^/[10-i2GeV-i2l 

D ' L J 


0.31(0.57) 


0.14(0.22) 


0.15(0.26) 


0.19(0.31) 


0.19(0.33) 


0.14(0.22) 


P8*^/[10-i6 GeV"^^] 


0.31(0.70) 


0.14(0.25) 


0.15(0.29) 


0.18(0.35) 


0.19(0.38) 


0.14(0.25) 


P*oV[10-'° GeV-20] 


0.34(0.89) 


0.14(0.28) 


0.16(0.35) 


0.19(0.42) 


0.20(0.46) 


0.14(0.30) 




0.63(3.75) 


0.24(0.80) 


0.27(1.12) 


0.33(1.36) 


0.34(1.54) 


0.24(0.92) 




Aihard = 5 GeV 
yUfac = 5 GeV 


/xsoft = 2.5 GeV 
/Xfac = 5 GeV 


/Xsoft = 2.5 GeV 
Pthard = 5 GeV 



Table 3: The theoretical moments P^^ for n = 4, 6, 8, 10, 20 and fixed as{Mz) = 0.118 and M^^ = 
4.8 GeV (Mf = 4.8 GeV) for various choices of the renormalization scales /isoft j /^hard 

and /Ufac- The 

two-loop running has been employed for the strong coupling. The values for P^^ in the pole mass 
scheme are slightly different from the numbers shown in Ref. [^ because of the correction in 02- 



We emphasise that the expression in Eq. (|T^) must be expanded consistently up to NNLO in the 
non-relativistic expansion together with the LO, NLO and NNLO contributions of Qn^i cind with Qn^2 
in order to achieve a proper cancellation of the correlations and the large NNLO corrections which are 
present in the pole mass scheme. Like in Ref. [Q] we keep the short-distance coefficient Ci in factorized 
form. 

Comparing Eqs. (|l^) and ( p!6|) we can easily see that the 15" scheme represents the most natural 
scheme one can use to avoid correlations and large corrections because it reduces the contribution 
from the lowest lying bottom-antibottom quark resonance, which represents the most important part 
of theoretical moments at large values of n. 

Before we turn to the statistical analysis it is quite interesting to study the impact of the 
transition to the IS scheme on the theoretical moments. In Tab. ^ we have displayed the values of 
P^ in the IS (pole) scheme for n = 4, 6, 8, 10, 20 and for various values of M^^ {Ml"^") and a,{M^) 
while the renormalization scales are fixed to //soft = 2.5 GeV and //hard = /^fac = 5 GeV. The numbers 
in Tab. ^ show that the dependence of the moments on the mass is quite strong in the IS* as well 
as in the pole mass scheme. However, in the 15 scheme the variation with respect to changes in the 
strong coupling is much weaker, in particular for larger values of n. To illustrate this we have also 
displayed the moments for n = 20, although this value it too high for the practical application. From 
this we can expect that the extracted values for are less strongly correlated to the strong coupling 
than the pole mass values obtained in Ref. [Q. This feature, however, will also make an independent 
determination of a precise value for the strong coupling from the T sum rules impossible in the IS" 
scheme, as we show in the next section. 

In Tab. ^ the theoretical moments are displayed for n = 4, 6, 8, 10, 20 for different choices of 
A*soft, Aihard and /x^c for as{Mz) = 0.118 and = 4.8 GeV (Mf"''' = 4.8 GeV) in the 15 (pole) 
mass scheme. As expected from the weak dependence of the moments in the 15 scheme on the strong 
coupling, we also find that the dependence of the moments on the soft scale is weaker in the 15 scheme. 
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Moment 


Ml^{Ml°'^)/[GeV] 


as{Mz) 




4.6 


4.8 


5.0 


0.11 


0.12 


0.13 


pth 

4,LO 


0.35(0.40) 


0.25(0.29) 


0.18(0.21) 


0.21(0.23) 


0.26(0.31) 


0.35(0.43) 


pth 

-'^4,NLO 


0.19(0.26) 


0.13(0.18) 


0.09(0.13) 


0.13(0.16) 


0.13(0.18) 


0.12(0.23) 


pth 

^4,NNLO 


0.27(0.39) 


0.19(0.28) 


0.14(0.20) 


0.17(0.22) 


0.20(0.30) 


0.26(0.46) 


TDth 


0.45(0.62) 


0.19(0.27) 


0.09(0.12) 


0.15(0.19) 


0.21(0.29) 


0.30(0.53) 


pth 

-^lO.NLO 


0.26(0.53) 


0.11(0.22) 


0.05(0.10) 


0.10(0.16) 


0.11(0.24) 


0.12(0.43) 


pth 

10,NNLO 


0.41(0.90) 


0.18(0.39) 


0.08(0.17) 


0.14(0.25) 


0.19(0.44) 


0.27(0.94) 




as{Mz) = 0.118 


Ml^ = 4.8 GeV 
(Mf = 4.8 GeV) 




/isoft = 2.5 GeV , /Xhard = Mfac = 5 GeV 



Table 4: The theoretical moments Pl^ and Pff^ at LO, NLO and NNLO for fixed /isoft = 2.5 GeV and 
/^hard = A*fac = 5 GeV for various values of (M^^"'") and a,(Mz) in the IS (pole) mass scheme. 
The two-loop running has been employed for the strong coupling. The values of and P^q given 
in units of 10"^ GeV~® and 10"^° GeV"^°, respectively. 



n 


4 


6 


8 


10 


P^^/[10~2n GeV-2"] 


0.22 


0.21 


0.23 


0.26 


(pex)T{15)/[io-2" GeV-2«] 


0.11 


0.13 


0.17 


0.21 


(pex)T{15)/pex 


0.48 


0.62 


0.72 


0.79 



Table 5: The experimental moments, -P^^, the contribution of the T(15') to the experimental moments, 
(^pex^T{i5)^ and the ratio (P^^)'^(i^) /P^^ for n = 4, 6, 8, 10. 



in particular for larger values of n. In the IS" scheme there is also some stabilization with respect to 
variations of the factorisation scale. The dependence of the moments on the hard scale, on the other 
hand, is comparable in both schemes. This is expected because in the transition from the pole to the 
IS* scheme only non-relativistic contributions are modified. Because the variation of the moments with 
respect to the renormalization scales is used as an instrument to estimate theoretical uncertainties in 
the mass extraction we can expect smaller uncertainties in the IS mass extraction compared to the 
results shown in Ref. 

In Tab. ^ the theoretical moment P^ and P^q are displayed at LO, NLO and NNLO in the 
IS (pole) scheme. Whereas in the pole mass scheme the NNLO corrections are all of the same size or 
even larger than the NLO corrections, they are significantly smaller in the IS scheme. This illustrates 
the improvement of the perturbative behaviour of the theoretical moments if the IS mass scheme is 
employed. We emphasise, however, that even in the IS scheme the NNLO corrections are not per se 
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small, particularly for large values of n. We believe that this is not a point of concern, because of the 
extreme dependence of the moments on the bottom quark mass, particularly for large values of n. It 
is the behavior of the results for the mass extracted from the theoretical moments which should be 
taken as the measure to judge the quality of the perturbative expansion. 

Finally, in Tab. ^ we have displayed for n = 4, 6, 8, 10 the values of the experimental moments 
P^^, based on Eq. (^) and the data given in Tab. |l|, and the contribution to P^^ coming from the 
T(15') in order to demonstrate the relative weight of the lowest lying resonance compared to the other 
resonances and the continuum. The contribution of the T(15) is between 48% (n = 4) and 79% 
(n = 10). This shows that the higher T excitations and the continuum constitute a significant part 
of the experimental moments. Thus, the fact that the value of M^^ which we determine in this work 
(see Eq. (22)) is compatible with is non-trivial. 



4 Numerical Results 

To obtain numerical results for the 15 mass we use the statistical procedure described in Ref. [§], 
which is based on the function 

(M^^, a.(Mz), //soft, /^hard, /^fac, a,or) = E {^n - ) (5"^)^™ ( P'rH - P^) ■ (19) 

{n},{m} 

{n} represents the set of n's for which the fit is carried out and is the inverse covariance ma- 
trix describing the experimental errors and the correlation between the experimental moments. The 
covariance matrix contains the errors in the T electronic decay widths, the electromagnetic coupling 
Oiem-, and the continuum cross section Vcont- The small errors in the T masses are neglected. The 
correlations between the individual measurements of the electronic decay widths are estimated to be 
equal to the product of the respective systematic errors given in Tab. || times the constant acor- In 
order to estimate the theoretical uncertainties in the mass extraction the renormalization scales /tgoft) 
A'hard, /"faci and the constant acor are varied randomly in the ranges 

1.5 GeV < /xsoft < 3.5 GeV 

2.5 GeV < /ihard < 10 GeV 

2.5 GeV < /xfac < 10 GeV 

< acor < 1 , (20) 



and the sets of n's 



{n} = {4, 5, 6, 7} , {7, 8, 9, 10} , {4, 6, 8, 10} (21) 



are employed. For each choice of the parameters /isoft; /"hard; A'fao and acor and each set of n's the 
value of M^^ is obtained for which is minimal. As in Ref. [Q] we carry out two types of fits. First, 
we consider a fit where M^^ and as{Mz) are determined simultaneously ("unconstrained fit"), and, 
second, we fit for M^'^ taking as{Mz) as an input ("constrained fit"). 

The result for the allowed region for M^^ and as{Mz) for the unconstrained fit using the full 
NNLO expressions for the theoretical moments is displayed in Fig. |^a. To illustrate the improvement 
coming from the NNLO corrections we have also displayed the result for the NLO theoretical moments 
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Figure 1: Result for the allowed region in the -Os plane for the unconstrained fit based on the 
theoretical moments at NNLO (a) and NLO (b). The dots represent points of minimal foi" ^ large 
number of random choices within the ranges ( |20| ) and the sets (^). Experimental errors are not 
displayed. The two-loop running has been employed for the strong coupling. 



in Fig. nib. The dots represent the points of minimal x ^ov a large number of random choices within 



the ranges ( poD and the sets (21). Comparing to the corresponding results in the pole mass scheme (see 
Fig. 10 and 12 in Ref. we find that the range of the mass values covered by the dots is smaller in the 
IS* scheme. For the NNLO moments the covered range for the IS mass is 4.67 GeV < M^'^ < 4.72 GeV 
(versus 4.76 GeV < M^°'^ < 4.85 GeV in the pole scheme [0). However, the range of Og covered 
in the 15" scheme is as large (slightly larger) as in the pole scheme at NNLO (NLO). This is not 
unexpected because of the reduced correlation of the theoretical moments to the strong coupling 
mentioned before. Obviously, the strong coupling determined from the T sum rules in the 15 scheme 
contains uncertainties which are much larger than the current world averages. We conclude that the 
T sum rules are not very competitive tool to determine the strong coupling. We therefore abandon 
the unconstrained fit also for the mass determination and turn to the constrained fit. 
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Figure 2: Result for the allowed M^^ values for a given value of at NNLO (a) and NLO (b). The 
dots represent points of minimal \^ for a large number of random choices within the ranges (pO[) and 



the sets (21), and randomly chosen values for the strong coupling 0^(2.5 GeV). Experimental errors 
at the 95% CL level are displayed as vertical lines. It is illustrated how the allowed range for M^^ is 
obtained if 0.114 < as{Mz) < 0.122 is taken as an input. The two-loop running has been employed 
for the strong coupling. 



In Fig. |2|a the allowed range for M^^ is displayed as a function of based on the NNLO 
theoretical moments. Fig. |2|b shows the result of the same analysis using the NLO theoretical moments 
only. Each dot represents the 15" mass for which the function is minimal for a given input value of a^, 
and for a random choice within the ranges ( pO|) and the sets (^). The statistical errors corresponding 
to 95% CL are below 15 MeV for all the dots displayed in Figs. |2|a and b. In contrast to the constrained 
fit results in the pole scheme (see Fig. 13 in Ref. Q) we find that in the IS scheme the correlation 
between the mass and the strong coupling is reduced considerably. Compared to the pole mass results 
the uncertainties are smaller by a factor three for the NLO as well as the NNLO analyses. Comparing 
Fig. ^ and b we also see that the NNLO contributions to the moments reduce the uncertainties of 
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the NLO analysis roughly by a factor two. Taking as{Mz) = 0.118 it 0.004 as input we arrive at 

= 4.71 ± 0.03 GeV (22) 

for the IS mass obtained from the NNLO theoretical moments, where also the experimental errors 
have been taken into account. This result is consistent with the unconstrained fit shown in Fig. |l|a. 
The NLO analysis yields M^^ = 4.70 it 0.06 GeV. Due to the small correlation of the IS" mass to 
the strong coupling the result does not change significantly if a slightly different choice of the central 
value and the error of as{Mz) is assumed. Eq. (^) represents the main result of this paper. We note 
that the choice of the range of the renormalization scale //soft has the main impact on the uncertainty 
in the IS mass determination. We have checked that allowing for larger values of //soft does not 
affect the results displayed in Figs. |2|. However, the uncertainties are larger if values of /Usoft smaller 
than 1.5 GeV are admitted. If a lower bound for /igoft of 1 GeV is chosen, the result for the 15 
mass reads M^^^ = 4.73 ± 0.05 GeV using the NNLO theoretical moments. However, we believe that 
//soft > 1.5 GeV is already a conservative choice. 

To conclude this section we would like to note that our method to extract the IS bottom 
quark mass based on the function for several moments (see Eq.(|l9|)) leads to smaller uncertainties 
than if the fit would be carried out for individual moments independently. This can be easily seen if 
the scale dependence of the moments shown in Tab. ^ is compared to their mass dependence shown 
in Tab. |2[ The function has a smaller scale dependence than the individual moments because it 
puts higher weight on linear combinations of the moments which are more sensitive to the relative 
size than to the absolute size of the moments. These linear combinations are determined from the 
entries of the covariance matrix (see Eqs. (81) and (82) of Ref. @) which account for the correlation of 
the experimental moments in Eq. (|l3| ) coming from their dependence on the T masses and electronic 
widths. 



5 Determination of the MS Mass and the Upsilon Expansion 

The IS mass is, by construction, optimally adapted to reduce the dependence of the theoretical 
moments on the strong coupling and the renormalization scale //soft; which govern the non-relativistic 
dynamics of the bottom-antibottom quark pair. By the same token, the 15 mass does not know much 
about large momenta above the inverse Bohr radius ~ M^a^. Thus, although the 15 mass can, like 
the MS mass, serve as a proper short-distance mass definition in its own right, it can only serve as 
a specialised short-distance mass definition for systems where the characteristic momentum scale is 
comparable to the inverse Bohr radius.^ For high energy processes, where momenta of order Mf, and 
higher are the characteristic scales, the MS mass is much better adapted. Thus, it is mandatory to 
determine the MS mass from the value of the 15 mass given in Eq. ([2^). We emphasise again that 
using specialised short-distance masses can reduce perturbative uncertainties because correlations are 

* By considering the fictitious case of inclusive semileptonic B decays with the radiation of d additional massless 
colour- singlet scalars at the weak vertex, where d is large, it was pointed out in Ref. that the characteristic scale 
for the inclusive B decay rates is Mi,/(5 + 2d) rather than M^. Because Mi,/{5 + 2d) ~ Adt as for d = 0, this might be 
an explanation for the drastic improvement of the convergence behaviour in the low order perturbative series for the 
inclusive B decay rates in the 15 scheme which has been reported in Ref. [0, 0. 
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minimised. This is a consequence of the fact that usually only a few low orders of perturbation theory 
are known. 

Great care has to be taken if the MS mass is expressed in terms of the IS" mass in order to 
ensure the proper cancellation of the most IR sensitive terms which exist in the perturbative series 
which expresses both masses in terms of M^"'*^. It has been shown in Ref. |11, ^ that the consistent 
way to relate the 15 and the MS mass to each other is the upsilon expansion. The upsilon expansion 
states that the LO, NLO and NNLO non-relativistic contributions in Eq. (|T|) are of order e, and 
e^, respectively, of the auxiliary parameter e = 1 mentioned already in Sec. |2|. On the other hand, 
in the relation which expresses the pole mass in terms of the MS mass the terms of order e, and 
in the upsilon expansion correspond to the one-, two- and three-loop contributions, respectively. 
When the MS mass is expressed in terms of the 15 mass we then have to use the expansion in e 
instead of a^. This means that at each order in e different orders of Og are mixed together. This 
unusual prescription can be understood from the fact that the most IR sensitive contributions in 
Eq. (H) (i.e. those contributions which involve the highest power of /3o in each order) involve powers of 
the logarithmic term L = ln(/x/Ci? Af^"'^). These logarithmic terms exponentiate at larger orders, 
X]i=o-^V^' ~ exp(L) = fi/Cp a-s M^"''^ , and effectively cancel one power of Ug |11, 12|. We will show 



below that this exponentiation is already very effective at NNLO in Eq. ([T|) making the use the upsilon 
expansion mandatory to achieve a reliable determination of the MS mass. 

The upsilon expansion states that we need the three-loop contribution in the relation between 
the MS and the pole mass in order to be able to make use of the known NNLO contribution ^^^^o 
m Eq. (I). For the purpose of illustration we will use in the following the large-/3o approximation 
for the yet unknown three-loop contributions. The relation between the pole and the MS mass then 



reads [23 



1 + 0.424 a, (mb) e + 0.158 a^(m;,) {po - 2.384) + 0.047 a^(m6)/3^ 



.(23) 



Combining Eq. (|^) with Eq. ( p3[) using the upsilon expansion up to we arrive at the following result 
for the MS bottom quark mass 



mb{mb) 



4.71 - 0.40 e - 0.11 - 0.04 ± 0.03 (SMf^] ± x 0.01 { 6a, 



GeV (24) 



where we have taken Eq. (22) as input for the IS mass and assumed as{Mz) = 0.118 it x 0.001 for 



the strong coupling. If the actual three-loop contribution in Eq. ( ^31 ) is just 10% smaller than the 
large-/3o approximation, the term in Eq. ( p^ ) is smaller by a factor of two. This illustrates the 
effectiveness of the exponentiation of the logarithmic terms mentioned above already at NNLO in 
Eq. (^). We emphasise that, because of this large cancellation, the MS mass would have a large 
(positive) systematic error if the NNLO contribution in Eq. (|l|) would be used without also including 
the three- loop contributions in Eq. (|2^). This systematic shift would amount to about -1-200 MeVj^ 
Thus at present stage we have to leave out all terms completely, and, assuming an uncertainty in 
as{Mz) of 0.004 {x = 4), we arrive at 

rnb(mfe) = 4.20 ± 0.06 GeV (25) 



This systematic error is included in results given in Ref. 
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for the MS bottom quark mass. To obtain the error in Eq. ( p5[ ) we have added all errors quadratically 
assuming a perturbative uncertainty due to ignorance of the exact size of the terms of 40 MeV. We 
believe that our optimistic assumption about the perturbative uncertainty is justified in view of the 
good quality of the large- /3o approximation which can be expected at the three-loop level in Eq. (23). 
In order to reduce the uncertainty in the MS bottom quark mass the three-loop contributions in 
the relation between the MS and the pole mass must be determined and the precision of the strong 
coupling has to be improved. 



6 The Low Scale Running and the Potential Subtracted Bottom 
Quark Masses 

In recent literature there have been two other proposals for specialised short-distance mass definitions 
which can be used for a stable and correlation-reduced bottom quark mass determination from the T 
sum rules. 

In Refs. |2^, the "low scale running mass' was proposed to subtract the low mo- 

mentum behaviour of the bottom quark self energy in the pole mass scheme.^ Like the 15 mass the 
"low scale running mass" was originally devised to improve the convergence of the perturbative contri- 
butions in inclusive B meson decays. Apart from the renormalization scale which governs the strong 
coupling the low scale running mass depends, in addition, on the cutoff uls which limits the momenta 
subtracted from the self energy. If this cutoff is adjusted to be close to the inverse Bohr radius the 
low scale running mass acts very similar as the IS mass. The relation between M^^ and the pole 
mass is known to two- loop order [26|, which corresponds to order in the upsilon expansion. In order 



to allow for compatibility checks for sum rule analyses and other bottom quark mass determinations 
within the framework of the low scale running mass scheme, it is useful to determine M^^ from the IS 
mass obtained in this work. Taking the large-/3o approximation for the yet unknown three-loop (order 
in the upsilon expansion) contribution, following Ref. ||5|, and using again the upsilon expansion 
one can determine M^^{ij,ls) from our result for the IS mass in Eq. (|2^). For ij-ls = ^ GeV the result 
reads 

Mt^{fiLs = I GeY) = 



4.71 - 0.09 e - 0.05 - 0.02 



0.03 (JM^J^) ± X 0.003 (ja^) ± 0.01 (/x) 



GeV , (26) 



assuming again as{Mz) = 0.118 it x 0.001 for the strong coupling as input and choosing /i = 5 GeV 
as the renormalization scale of the strong coupling. The third uncertainty in Eq. (|2^) arises from 
the variation of this renormalization scale of the strong coupling between 2.5 and 10 GeV. In the 
relation between the IS and the low scale running mass it is reasonable to choose the renormalization 
scale of order the bottom quark mass because the physical effects involving the inverse Bohr radius 
are eliminated. [Allowing for /i = 1.5 GeV would not change the result shown in Eq. (p6|).] The 
uncertainty from in Eq. (|26|) (and also in Eq. (|28|)) is smaller than in the corresponding uncertainty 

^ In this language the IS mass just subtracts half the perturbative binding energy of a ^Si bottom-antibottom quark 
bound state in the pole mass scheme. 
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in Eq. (|24| ) because the order e term is much smaher in Eq. (|2q ) (and in Eq. (psD). The relation 
between M^^ and is very well behaved illustrating the fact the both are short-distance masses. 
The size of the scale uncertainty is compatible with the size of the term. Assuming an uncertainty 
of 0.004 (x = 4) for as{Mz) and taking into account that the exact e'^ contribution is not yet known 
we arrive at 

Mfe^^(/UL5 = 1 GeV) = 4.57 ±0.04 GeV (27) 

for the low scale running mass at 1 GeV. As in the previous section we have combined the error 
coming from M^^ , a^, and a perturbative uncertainty of 20 MeV quadratically. The result for the 
low scale running mass in Eq. (27) is compatible with a recent determination of M^^{^ls = 1 GeV) 



which was also based on the T sum rules at NNLO (M^^-^ = 4.56 ± 0.06 GeV) g. The analysis 
of Ref. m, however, is different from ours with respect to several aspects. In Ref. the energy 
denominators in the cross section R{s) have been resummed for energies close to the resonances, and 
the final dispersion integration (^) has been carried out numerically, which means that the moments 
calculated in Ref. contain certain non-relativistic contributions beyond NNLO. In this work the 
energy denominators are not resummed and all integrations are carried out analytically, consistently 
dropping all higher order non-relativistic contributions beyond NNLO (see Ref. Q for details). This 
means that our analysis only uses global duality arguments, whereas Ref. also relies on the validity 
of local duality. In addition, in Ref. Q larger values of n (n = 14 — 18) and larger values of the soft 
scale /isoft (A*soft = 2 — 4.5 GeV) where used, and M^^ was determined from extracting values for the 
pole mass first and converting them into numbers for the low scale running mass afterwards. In our 
analysis the pole mass has been eliminated completely. We also point out that in Ref. |^] the values 
of M^^ were obtained by fitting individual moments whereas in this work a x function based on 
several moments was employed (see the comment at the end of Sec. §). The compatibility of Eq. (27) 



with the result obtained in Ref. might serve as an argument that the methodical differences in this 
analysis and the one in Ref. |^] do not affect the mass determination. 

In Ref. the "potential subtracted" ,Mft^^, mass was proposed. It subtracts the low mo- 



mentum contribution of the static potential in the pole mass scheme |15| , |l6| , 17j. Like the low scale 
running mass the potential subtracted mass depends on a cutoff, fips- For ^ps = ^fJ-LS the low scale 
running and the potential subtracted mass are approximately equal. This equivalence is based on 
the universality of the most infrared sensitive contributions contained in the corresponding subtrac- 
tions. So far there has not been any sum rule analysis which has attempted to extract the potential 
subtracted bottom quark mass. However, it is useful to determine M^"^ from the 15 mass value 
obtained in this work to allow for cross checks with possible future work on this subject. Because the 
potential subtracted mass is based on the static perturbative potential, which is known to NNLO in 
the non-relativistic expansion (i.e. in the upsilon expansion), the perturbative uncertainties in the 
relation between potential subtracted and 15" mass are smaller. Starting from our result for the IS 
mass, Eq. (|2^), the potential subtracted mass at the subtraction scale hps = 2 GeV reads 



M^if^PS = 2 GeV) 



4.71 - 0.13 e - 0.04 - 0.01 



Om(^6Mt^^ ± X 0.003 (fcs) ± 0.01 ( 



GeV , (28) 



where we have again assumed as{Mz) = O.llSitxO.OOl as an input for the strong coupling. The third 
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uncertainty arises from the variation of the renormalization scale in the strong couphng between 2.5 
and 10 GeV. The relation between M^^ and is very well behaved illustrating the fact that M^^ 
is also a short-distance mass. Assuming again an uncertainty of 0.004 (x = 4) for ^^(M^), taking a 
perturbative uncertainty of 10 MeV and adding all errors quadratically we arrive at 

mP{iips = 2 GeV) = 4.53 ± 0.03 GeV (29) 

for the potential subtracted mass at 2 GeV. 

To conclude this section we note that it is possible to tunc the cutoff scales in the low scale 
running and the potential subtracted masses, at each order of perturbation theory, such that the 
correlations of the mass to other parameters are minimal. Such a fine-tuning approach, however, is 
illegal because from the conceptual point of view any cutoff scale close to the inverse Bohr radius is 
equally well suited. This "scale-ambiguity" is a specific characteristic of the low scale running and 
the potential subtracted mass and not inherent in the IS" mass, because the latter is defined through 
a bound state mass, which is a physical quantity. We therefore consider the IS* scheme as the most 
natural mass scheme to be used to describe processes involving heavy quark-antiquark pairs in the 
non-relativistic regime. We emphasise, however, that the L5, PS and IS mass improve the quality 
of the mass extraction from the T sum rules all by the same mechanism, and that the most infrared 
sensitive contributions contained in the subtractions are equivalent in all schemes. In the IS scheme 
the cutoff which is visible in the definition of the LS and the PS masses is provided in a natural way 
by the width of the wave-function of the bottom-antibottom quark bound state. No matter which 
scheme is used the final answers for the MS bottom quark mass should be compatible. 



7 Conclusions 

In this work we have extracted the IS* mass, M^^ , from sum rules which relate the masses and electronic 
decay widths of the T mesons to large-ra moments of the vacuum polarization function, which have 
been calculated at NNLO in the non-relativistic expansion. is defined as half the perturbative 

mass of a J^^ = 1 , bottom-antibottom quark bound state, i.e. it represents the perturbative 
contribution of half the mass of the T(IS') meson. However, the latter information is not used in this 
work and M^^^ is considered as a fictitious mass parameter which is determined solely from the sum 
rule analysis. The result for the IS" mass reads M^^ = 4.71 ± 0.03 GeV, if as{Mz) = 0.118 ± 0.004 is 
taken as an input for the strong coupling, and has only a small correlation to the choice of ag ■ The close 
proximity of this result to Mx(is)/2 is non-trivial because higher excitations and the hh continuum 
have a non-negligible contribution to the sum rules for the values of n used in this work. The fact 
that, within the errors, is equal to Mx(i5)/2 = 4.730 GeV indicates that the non-perturbative 
contributions in the T(IS') mass are probably small. Compared to an earlier work which has been 
carried out in the pole mass scheme by the same author we find an improved perturbative behaviour 
of the moments and a much smaller residual dependence on the renormalization scale which governs 
the non-relativistic dynamics of the bottom-antibottom quark pair. The IS bottom quark mass is a 
short-distance mass (i.e. it does not have an ambiguity of order Kqcd) and can be used as a mass 
parameter in its own right. It is well adapted to processes where the characteristic scale is below the 
bottom quark mass and has been successfully applied earlier to inclusive B mesons decays, where it 
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leads to a significant improvement in the behaviour of the perturbative series describing the inclusive 
widths. [0, m 

Because M^^ is a short-distance mass it has a well behaved perturbative relation to the MS 
mass, which is adapted to processes where the characteristic scale is of order the bottom quark mass 
or higher. In order to relate the IS* mass to the MS mass definition it is mandatory to use the upsilon 
expansion |11, 12 1 in order to ensure the cancellation of large infrared-sensitive contributions present 
in the perturbative series which relate both masses to the pole mass. Using the presently known 
two-loop contributions in the relation of the MS and the pole mass (i.e. using the upsilon expansion 
up to order e^) and assuming as{Mz) = 0.118 it 0.004 we arrive at mb(mf,) = 4.20 it 0.06 GeV for the 
MS mass. The error in the rnbijnb) can be reduced once the three- loop contributions in the relation 
between the MS and the pole mass are calculated and if the precision in as is increased. 
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